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ABSTRACT 

We address the issue of the production of Lyman continuum photons by T Tauri stars, in an 
attempt to provide constraints on theoretical models of disc photoionisation. By treating the 
accretion shock as a hotspot on the stellar surface we show that Lyman continuum photons 
are produced at a rate approximately three orders of magnitude lower than that produced by a 
corresponding black body, and that a strong Lyman continuum is only emitted for high mass 
accretion rates. When our models are extended to include a column of material accreting on 
to the hotspot we find that the accretion column is extremely optically thick to Lyman con- 
tinuum photons. Further, we find that radiative recombination of hydrogen atoms within the 
column is not an efficient means of producing photons with energies greater than 13.6eV, and 
find that an accretion column of any conceivable height suppresses the emission of Lyman 
continuum photons to a level below or comparable to that expected from the stellar photo- 
sphere. The photospheric Lyman continuum is itself much too weak to affect disc evolution 
significantly, and we find that the Lyman continuum emitted by an accretion shock is simi- 
larly unable to influence disc evolution significantly. This result has important consequences 
for models which use photoionisation as a mechanism to drive the dispersal of circumstellar 
discs, essentially proving that an additional source of Lyman continuum photons must exist if 
disc photoionisation is to be significant. 

Key words: accretion, accretion discs - circumstellar matter - planetary systems: protoplan- 
etary discs - stars: pre-main-sequence 



1 INTRODUCTION 

The evolution and eventual dispersal of discs around young stars 
is an important area of study, as it provides constraints on the- 
ories of both star and planet formation. It is now well estab- 
lished that the majority of young stars at ages of ~ 10 6 yr have 
circumstellar discs which are o ptically th ick at optical and in- 
frared wavelengt hs IStro m et alll98^ : lKenvon & Hartmanrll995t 
lHaisch. Lada & Ladal l200lF ~ These discs are relatively massive, 
with masses of the order of a few percent of a solar mass 
jBeckwith et alJll99(jlEisner & Carpented f2003). However by an 
age of ~ 10 7 yr most stars are no longer seen to have such 
ma ssive discs, a lthough low- mass "debris discs" may r emain 
(eg. lMannings & Sargenil997tlWvatt. Dent & Greavej2003l) . The 
mechanism by which these discs are dispersed remains an impor- 
tant unsolved question. 

One fact which is clear, however, is that a large fraction 
of the mass from the disc will eventually be accreted on to 
the central (proto)star. The inner edge of the accretion disc is 
typically truncated by the magnetosphere at a radius of a round 
5R, iMever. Calvet & Hillenbrand 1 19971: ICalvet & Gullbrind 
1998) and so material falling from the disc on to the central star 
can attain an extremely high velocity, resulting in a so-called 



email: rda@ast.cam.ac.uk 



"accretion shock" when this material impacts upon the stellar sur 
face. Existing models of the accretion sho c k j Calve^^^ullbrin 
1 19981 IGullbring etaflll998t lLamzij[T99sl IGullbring et al]|200l 
have paid a great deal of attention to the emission in the ul- 
traviolet (1000-3000A) and visible (3500-7000A) wavebands, 
comparing theoretical predictions to observed spectra. Extremely 
good models have been constructed, and these emission spectra 
are now well understood. However, very little attention has 
been paid to the emission shortward of the Lyman break (< 
912A), primarily because absorption by interstellar Hi makes it 
impossible to observe young star s in this wavelength re gime. 
Recen t the oretical studies Ici arke, Gendrin & Sotomavor 

[2001; (Matsuvama, Johnston e & Hartmannl 2003; 
Armitage. Clarke & Pallal 12003) have suggested that photoioni- 
sation by the central object may play an important role in disc 
dispersal, and so the emission shortward of the Lyman break has 
become important. 

Currently the origin of photoionising radiation from young 
stars such as T Tauri stars is unclear, and even the m agni- 
tude of such emission is p oorly constrained 1 Gah m*et alJll979t 
Imhof f & AppenzelleJ 19871) . To date, models of disc photoionisa- 
tion have used either a co nstant ionising flux (eg. Hol lenbach et alJ 
1 1994 IClarke et al J 1200 ll) . assumed to be chromospheric in ori- 
gin, or modelled the accretion-driven flux simply as a con- 
stant temperature hotspot on the stellar surface, emitting as a 
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blackbody iMats uvama et alJl2003l) . The latter produces an ion- 
ising flux which is proportional to the mass accretion rate, 
thus decreasing dramatically with ti me as the ac c retion rate 
falls, and so naturally the models of Iciarke et alJ 1200 ll) and 
IMatsuvama et al J 120031) h ave produced markedly diff erent results . 
These models iHollenbach et alJ Il994t iRichling & Yorkd 11997$ 
iHollenbach. Yorke & Johnstonell200(]|) also find that the mass-loss 
rate scales approximately with the square root of the ionising flux, 
and so very high ionising fluxes are required to influence disc evo- 
lution significantly: an ionising flux of 10 41 photons s -1 will drive 
mass-loss from the disc at around the 10~ 10 A/eyr -1 level. 

In this paper we study the issue of the ionising flux gener- 
ated by an accret ion shock. The simplifie d model of the accretion 
shock adopted by Mat suvama et a l. 1 2003) neglects two key points, 
which we address in turn. Firstly, it seems likely that the Lyman 
continuum emission from such a hotspot would resemble a stel- 
lar atmosphere rather than a blackbody; whilst the two are almost 
identical at longer wavelengths, photo-absorption by Hi provides 
a strong suppression of the flux shortward of 912A. Secondly, the 
photons emitted by the accretion shock must pass through the col- 
umn of material accreting on to it in order to interact with material 
in the disc, and again we would expect photo-absorption by Hi in 
the column to supress the Lyman continuum significantly. In or- 
der to address these issues we have modelled this process in some 
detail. In section|2|we investigate the effect of replacing the black- 
body hotspot with a more realistic stellar atmosphere. In section|3| 
we investigate the effect of passing these photons, from both the 
blackbody and the stellar atmosphere, through a column of accret- 
ing material. In section[4]we discuss our results and the limitations 
of our analysis, and in section|5|we summarise our conclusions. 



2 STELLAR ATMOSPHERES 

The photoionisation models of IMatsuvama et al] 120031) model the 
ionising photons as follows. They assume that the flux from the 
accretion shock can be modelled as a constant temperature hotspot, 
and adopt a blackbody spectral energy distribution at a temperature 
of T = 15, 000K. They assume that half the accretion luminosity 
is radiated by this hotspot, and so for a star of mass A/» and radius 
_R* the accretion shock luminosity L is given by: 



GM,M d 
2R„ 



AosbT 



(1) 



where A/d is the rate of mass accretion from the disc, A is the area 
of the hotspot, and <jsb is the Stefan-Boltzmann constant. As the 
temperature is constant the rate of ionising phot ons <ft a is s imply 
propor tional to the accretion rate Afd • In addition, Mat suvama et alJ 
(2003) add a further contribution to the ionising flux from the stel- 
lar photosphere, neglecting the poorly-constrained chromospheric 
contribution, at a constant rate of <I?p = 1.29 x 10 31 photons s _1 , 
with the total ionising photon rate given by <!> a + <!>p. Our assertion 
is that, given the high density of atomic hydrogen, it is extremely 
unlikely that such a hotspot would radiate as a blackbody. A spec- 
trum akin to a stellar atmosphere, showing a significant "Lyman 
edge", seems much more likely. As a result, we first consider the 
strength of this effect. 



2.1 Constant temperature 

Our first models simply involve substituting model stellar atmo- 
spheres in place of the blackbody emission in equation^ We have 
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Figure 1. Spectra of light incident on and emitted by a constant density ac- 
cretion column with a covering factor of 5% and a height of IRq . Note the 
precipitous drop in the emitted spectra at the Lyman break, due to photoab- 
sorption by Hi. The deep absorption feature at 1215A is Lya. 



adopted the same stellar parameters as [Matsuvama et alJ 120031) 
(R — IRq, M — IMq), and similarly adopted a constant hotspot 
temperature of 15,000K. The luminosity L scales with Md in the 
same manner as in equ ation [J] We have utilised the Kurucz model 
atmospheres jKuruc3[l992) (which have been incorporated into 
the CLOUDY code) for this temperature and surface gravity. The 
model atmospheres do not deviate significantly from the blackbody 
at longer wavelengths, but are some 3 orders of magnitude less 
luminous than the corresponding blackbody at wavelengths short- 
ward of the Lyman limit at 912A, due to absorption by Hi. Fig.Q 
compares the blackbody and Kurucz spectra. (The apparent lack of 
emission lines in FigQis an artefact of the relatively large bin-width 
used within the code. The use of such large wavelength bins results 
in line fluxes which are negligible in comparison to the continuum.) 
FiglUplots the ionising fluxes as a function of accretion rate, and 
we see that the stellar atmosphere hotspot produces ionising pho- 
tons at a rate that is a factor of 1100 less than that obtained from 
the blackbody model. 

2.2 Constant area 

Another consider ation is that of the hotspot area. The 
IMatsuvama et alJ 120031) blackbody formulation described in 
equation Q uses a hotspot temperature which remains constant 
for different mass accretion rates, implying a hotspot size which 
decreases as the accretion rate drops. The accreting material is 
thought to be channelled on to the magnetic poles as it falls on 
to the stellar surface, and unless the topology of the magnetic 
field varies systematically wth the accretion rate, the hotspot area 
should remain approximately constant with time. If the hotspot 
area A is constant then we would expect, from equation Q the 
hotspot temperature to vary as T a A/ d The total luminosity 
of a model atmosphere is very similar to that of a blackbody, 
and so we adopt this relationship for the stellar atmospheres also. 
We re-evaluated the model atmospheres described in section |2~T1 
keeping the scaling luminosity proportional to Afd, but now with 
the temperature given by: 



T = 



GM, 



2R*Aa SB 



Mr 



1/4 



(2) 
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Figure 2. Ionising flux rates from different accretion shock models. The 

rate from a simple model atmosphere is some 3 orders of magnitude less 

■ 1 /4 

than that from a blackbody, and the constant area case, with T oc M, 
decays more precipitously still. 



The result of this is shown in Fig|2| the drop-off in the ionising 
photon rate is much more precipitous than in the constant tem- 
perature case, with only very high mass accretion rates, greater 
than 10 -7 Moyr -1 , producing ionising photons at greater than the 
photospheric rate. In fact, the relationship between the hotspot 
area and the mass accre tion rate is not well understood, and 
ICalvet & GullbrinsU l998) even found observational evidence for a 
hotspot area which increases with Md (which would imply an even 
steeper decline in the ionising flux). However, we go on to show 
that the presence of an accretion column above the hotspot is by far 
the dominant factor in controlling the ionising photon rate, and so 
the exact details of the hotspot area are not of great significance. 



3 ACCRETION COLUMNS 

The other issue which affects the emitted ionising flux is the as- 
sumed presence of a column of accreting material directly above 
the hotspot. This material will absorb Lyman continuum photons 
through photoionisation of Hi, and so a large attenuation of the ion- 
ising flux is exp ected. Adopting the photoionisation cross-section 
from lCoxl 120001) of CT13.60V = 6.3x 10 _18 cm 2 , indicates that any 
column density greater than 5xl0 18 cm -2 will result in an atten- 
uation of the incident flux by a factor of > 10 13 , enough to re- 
duce any incident ionising photon rate to below photospheric lev- 
els. The density of the inf a iling material is of order 5x 10 12 cm -3 
ICalvet & Gullbrina [l998t) and so this results in an attenuation 
length of order 10 6 cm (1O _5 7?0). As a result, the only Lyman con- 
tinuum photons which can be emitted, at any significant rate, by an 
accretion column must be due to radiative recombination of hydro- 
gen in the column, the so-called diffuse continuum. 

In order to investigate this effect further we have constructed 
models of the accreti on column using the CLOUDY photoionisation 
code jFerland [l996). These models consist of a uniform central 
source, radiating either as a blackbody or a stellar atmosphere, and 
an accretion column. The central source emits in the radial direction 
only, and the accretion column covers a constant solid angle. Thus, 
by a simple linear subtraction of the emission not incident on the 
column we can treat the column as if it were illuminated solely by 
a hotspot at its base with a flux in the radial direction. In reality 



such a hotspot would produce some lateral component of flux near 
to its edges: this is discussed in section |4] The accretion column 
has a solar chemical composition (although as the dominant effect 
is photoionisation of Hi there is almost no dependence on chemical 
compositio n) and covers a small fra ction / of the stellar surface. 
Following lCalvet & Gullbrinl (l 998 ) we adopt a "free-fall" scaling 
density of the form: 



n H (0) = 5.2 x 10 12 cm" 3 



M d 



10- 8 M Q yr- 1 
/ R* 



\2R Q 



0.5M© ) 
3/2 / , 

aoi 



(3) 



The radial behaviour of the density is obtained by assuming 
that the material falls along magnetic field lines, in a man- 
ner consistent with stand ard magnetospheric accretion models 
(eg. iGhosh & Lamb l ll978h . Therefore at a given radius the prod- 
uct of the field strength and the column cross-sectional area is a 
constant. Assuming a dipole magnetic field we have B oc R~ 3 , 
and therefore the cross-sectional area is proportional to R 3 . For a 
column with cross-sectional area A' and mass density p, mass con- 
servation requires that: 



pA'vg — constant 



(4) 



The number density nn oc p, and the free-fall velocity vg oc 
ir 1/2 , and so we have a density scaling law of: 



riH oc R 



-5/2 



(5) 



with the condition in equation[3]used to fix the scaling constant. 

CLOUDY allows us to evaluate the continuum and line emis- 
sion from the top of the accretion column, which is a combination 
of the continuum incident on the bottom of the column, attenuated 
by the column, and the diffuse emission from the heated column. It 
does not allow direct evaluation of the emission from the "sides" of 
the accretion column, which may be significant and is discussed in 
section|4| 

Again, we have adopted R = IRq and M = IMq, and have 
constructed models of these accretion columns for a broad range 
of accretion rates, hotspot areas and column heights. Initially both 
the blackbody and stellar atmosphere hotspot formulations were 
used, with the "constant temperature" formalism used as it pro- 
vides the greatest ionising flux to the column and can be treated 
as a limiting case. However, as discussed above, the incident Ly- 
man continua from both are extinguished over a very short length 
scale, and so the only Lyman continuum emission which emerges 
from the columns is the so-called "diffuse" emission due to the ra- 
diative recombination of atomic hydrogen. As seen in Fig0 the 
Lyman continua emitted by the columns are identical in both cases, 
and depend only the nature of the column rather than the spectrum 
of the illuminating hotspot, as both hotspot spectra have the same 
bolometric luminosity. As a result, only the more realistic stellar 
atmosphere models were used for the remainder of the cases. 

3.1 Results 

The results of the simulations described above are presented in 
Figs0 & [4] Fig|3| shows the incident and transmitted spectra 
for a typical case at a variety of different column heights, and 
Fig|4] shows the dependence of the emergent ionising flux on col- 
umn height. CLOUDY is only valid for temperatures greater than 
~3000K - for temperatures lower than this the thermal solutions are 
no longer unique - and for large column heights or cases with little 
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Figure 3. Incident (solid line) and transmitted (dashed lines) spectra for 
model accretion column with = 10 — 9 MQyr _1 and a covering frac- 
tion of 1%. The spectra shown are for column heights of 0.2, 1.0, 2.0, 3.0 
and 4.0i?Q (brightest to faintest respectively). 

heating (ie. low accretion rates) the temperature dropped below this 
critical value. The models were not pursued beyond this point, with 
3000K used as the temperature limit of the calculations. As seen 
in Fig|4] the emergent ionising flux, as expected, decreases with 
both decreasing accretion rate and with increasing column height. 
Smaller hotspots result in higher density of material in the accretion 
column for a given accretion rate, and so tend to produce slightly 
larger ionising fluxes. However these variations are small relative 
to those due to variations in column height or accretion rate. 

The most important point, however, is that the emergent ion- 
ising photon rates for all of the columns are less than the photo- 
spheric value of 10 31 photons s _1 . Further, the photospheric value 
is itself some 10 orders of magnitude lower than the rate required 
to influence disc evolution significantly. This means that the accre- 
tion columns we have modelled cannot emit ionising photons at a 
rate that will be significant in disc photoionisation models, for any 
choice of parameters. 



4 DISCUSSION 

There are obvious caveats to these models. The first, and most 
significant, is that the CLOUDY code is only completely reliable 
at densities less than 10 13 cm -3 ; it is prone to numerical prob- 
lems at higher densities, and for cases of high Md and small / 
the density in our models can exceed this value. However, the 
main uncertainties at these densities are regarding the treatment 
of heavy elements. The dominant effects in the regime in which 
we are interested are photoionisation and recombination of atomic 
hydrogen, and these processes are treated reliably by the code. 
However in order to compensate for this numerical problem we 
were forced to limit the density artificially to have a maximum 
value of 5 x 10 13 cm -3 . This affected the 3 models with the high- 
est densities (M d = 1 x lO~ 7 M yr _1 , / = 0.01, 0.05 and 
M d = 5 x lO _8 M yr _1 , / = 0.01), and the reduction in density 
results in these models under-estimating the ionising flux some- 
what. This does introduce some uncertainty into our models, but 
we consider the effect to be small relative to the gross effects which 
dominate the calculations. 

A further caveat regards the issue of non-radial emission, both 
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Figure 4. Ionising photon rates emitted by accretion columns of different 
heights, for columns with covering factors (and hotspot areas) of: a) 1%; b) 
5%; c) 10% . From top to bottom in each plot the curves are for accretion 
rates of 1 x 10" 7 , 5 x 10~ 8 , 1 x 5 x 10~ 9 , 1 x 10~ 9 , 5 x 

10 — 10 & 1 x 10~ lo M0yr -1 respectively. The smallest column height 
in each case was taken to be half the hotspot radius. Where the curves stop 
abruptly it is because the temperature in the column fell below the limiting 
temperature of 3000K, and for the same reason the last curves (Md — 
1 x 10 — lo M0yr — 1 ) in plots b) and c) have been omitted completely. 
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from the hotspot and from the "sides" of the column, as mentioned 
in section [5] The ionising photon rates from our models are those 
emitted from the top of the accretion columns only, and neglect 
emission from the sides of the columns. Further, the incident flux 
from the hotspot is assumed to be purely radial. The Lyman contin- 
uum emission in which we are interested arises from the radiative 
recombination of hydrogen atoms, an intrinsically isotropic emis- 
sion process, and so emission from the sides of the column could 
be significant. However, as the recombination process is isotropic 
it is reasonable to assume that the Lyman continuum emitted from 
the sides of a column will be comparable to that emitted from the 
top of a column that is truncated at a height equal to its diameter. 
FigS shows that the diffuse Lyman continuum from the top of the 
column decreases dramatically as the column height increases (and 
the density decreases). Consequently the Lyman continuum emis- 
sion from the sides of the column will only be significant over a 
distance comparable to the hotspot diameter: the sides of the col- 
umn only produce a significant Lyman continuum near to the stellar 
surface. 

As noted in section[3] Lyman continuum photons incident on 
the columns will be attenuated by a factor of 10 13 over a distance 
of approximately 1O~ 5 7?0. As a result of this, the majority of the 
non-radial Lyman continuum photons emitted by an isotropically 
emitting hotspot will be absorbed by the column. The only such 
photons not absorbed will be those emitted within a fraction of an 
attenuation length of the edge of the hotspot, in a direction away 
from the centre of the hotspot. A hotspot covering 1% of the stellar 
surface has a radius of Q.2Rq, so less than 10~ 5 of the hotspot 
photons are emitted within lO _5 i?0 of the hotspot edge. Given this 
fact, and also the behaviour of the "raw" hotspot emission described 
in Fig[2| we neglect this effect. 

The net result of these two simplifications is that in reality 
the emission from the bottom of any accretion column will dom- 
inate the Lyman continuum emitted by the column, and accretion 
columns of any height will emit ionising photons at a rate compa- 
rable to that provided by the lower part of the column (the left-hand 
end of the curves in Fig|4). This will increase the largest calculated 
ionising photon rates by a small geometric factor but still cannot 
increase their flux to significantly greater than 10 31 photons s" 1 ; 
the photospheric emission will still dominate the overall Lyman 
continuum. More importantly, much higher ionising fluxes (of the 
order of 10 41 phot ons s -1 ) are required to have a significant effect 
on disc evolution lOarke et all200ltlMatsuvama et all2003t) . and 
the uncertainties caused by the approximations we have made are 
negligible compared to this 10 orders-of-magnitude difference. 

Further simplifications used in our models regard the geome- 
try of the accretion column. Our models use a column with a con- 
stant covering factor - essentially a truncated radial cone - and so 
the area of the outer surface at a given radius R is proportional to 
R 2 . However in reality the column is channelled by the magneto- 
sphere and, as explained in section [3] has an area proportional to 
R s ; our model accretion columns are somewhat less flared than we 
would expect to see in reality. However the difference between the 
two is only significant at large radii; small column heights, which 
provide the highest ionising fluxes, will not show significant devi- 
ation between the two cases. Again, the net result of this is that we 
probably under-estimate the ionising flux by a small factor, but not 
by enough to alter the results significantly. 

As our accretion columns are radial cones they do not bend 
to follow the magnetosphere, as expected in more rea listic mag- 
netospheric accretion models (eg. lGhosh & Lambll978l) . In such a 
model the column would bend over to meet the accretion disc, with 



a curvature dependent on both the latitude of the hotspot and the 
strength of the magnetosphere. However, as discussed above, the 
emission from the bottom part of the accretion column dominates 
over that from the upper parts (those affected by this curvature), so 
this simplification will not affect our results significantly. 

There is also the issue of an infall velocity, which our models 
do not address. In reality the accretion column will be falling to- 
wards the stellar surface at close to the free-fall velocity, which can 
be several hundred kms" 1 , and so this could modify the absorbing 
effect of the cloud. However the infall velocity is much less than 
the speed of light, and there are no strong emission lines near to the 
Lyman break, and so we consider the impact of this effect on the 
emitted Lyman continuum to be negligible. 

It should be noted that our model makes no predictions as to 
the behaviour of the photons emitted at wavelengths longward of 
the Lyman break. As shown in FigQ the emission from the top of 
the column longward of the Lyman break is essentially identical to 
that from the hotspot at the base of the column. We approximate 
the accretion shock crudely, and so are not able to fit our models to 
o bserved spectra in a manner similar tolCalvet & Gullbrind 1 19981) 
or ljohns-Krull. Valenti & Linskvl fcOOOh . However the Lyman con- 
tinuum emitted by our columns is insensitive to variations in the 
hotspot spectrum, due to the high optical depth of the columns to 
Lyman continuum photons. Consequently we find that any reason- 
able accretion-shock model will produce a similar Lyman contin- 
uum, and that this Lyman continuum is independent of the emission 
at longer wavelengths. 

We have adopted stellar parameters of R = lR p> and M = 
IMp) to pr ovide direct comparisons with the models of Clar ke et"al] 
1 200 1 ) and Matsuvama et al. 1 2003 ). However in the case of T Tauri 
sta rs a radius of 2Rp, a nd a mass of O.5A/0 would be more realis- 
tic iGullbring et al.ll99a) . The result of this will be that our models 
over-estimate the ionising flux somewhat, due to both reduction in 
the energy released by accretion on to the stellar surface and also 
due to a reduction in the star's surface gravity. Once again, how- 
ever, it is unlikely that these factors are significant in comparison 
to the gross effects we have already considered. Similarly, the use 
of the "constant temperature" hotspot to heat the column probably 
over-estimates both the ionising flux and heating provided by the 
hotspot, and thus over-estimates the diffuse Lyman continuum. In 
effect we have constructed a "best-case" model, designed to pro- 
duce the maximum ionising flux, and still found the ionising flux 
to be less than that emitted by the stellar photosphere. It seems ex- 
tremely unlikely that any conceivable accretion column could pro- 
duce ionising photons at a rate significantly greater than this. 



5 SUMMARY 

In an attempt to provide some constraints on the nature and mag- 
nitude of the ionising continuum emitted by T Tauri stars, we have 
constructed models which treat the accretion shock as a hotspot on 
the stellar surface beneath a column of accreting material. We have 
modelled these columns for a variety of different accretion rates, 
hotspot sizes and column heights, and have found that: 

• A hotspot radiating like a stellar atmosphere radiates ionising 
photons at a rate some 3 orders of magnitude less than the corre- 
sponding blackbody. 

• A constant area hotspot radiating like a stellar atmosphere can 
only emit ionising photons at greater than photospheric rates for 
mass accretion rates greater than lO _7 M0yr _1 . Such accretion 
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rates are nea r the upper limi t of the rates deriv ed from observations 
lHartmann et alJl998tljohns-Krull et alJ200Ci) . 

• Photoionisation of neutral hydrogen in the accretion column 
attenuates the Lyman continuum from any hotspot to zero over a 
very short length scale. The ionising photons which do emerge are 
due to radiative recombination of hydrogen atoms in the column, 
and the rate of ionising photon emission is less than the photo- 
spheric level for all of the accretion columns we have modelled. 

In short, we find that accretion shocks and columns are extremely 
unlikely to produce Lyman continuum photons at a rate signifi- 
cantly greater than that expected from the stellar photosphere. The 
photospheric level itself is some 10 orders of magnitude below the 
rates required for photoionisation to affect disc evolution signifi- 
cantly, and so it seems that the Lyman continuum emitted by an 
accretion shock will not be large enough to be significant in disc 
photoionisation models. These models have provided an attrac- 
tive explanation of some ob served disc properties (eg. Clar ke et"all 
2001; Armitase et al. 2003) but we have shown, as suggested by 
Clarke et al. 1 2001), that they must be powered by something other 
than the accretion-shock emission. In order to produce ionising 
photons at a rate large enough to enable photoionisation models to 
provide a realistic means of disc dispersal, the central objects must 
sustain a source of ionising photons that is not driven by accretion 
from the disc. 
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